TIFR/TH/00-66 



Electron and Neutron Electric Dipole Moments in the 
Focus Point Scenario of SUGRA Model 



Utpal Chattopadhyay (a) , Tarek Ibrahim^ and D.P. Roy (a) 

Department of Theoretical Physics 
Tata Institute of Fundamental Research, Mumbai - 400 005, India 

^ Department of Physics, Faculty of Science 
University of Alexandria, Alexandria, Egypt 



Abstract 

We estimate the electron and neutron electric dipole moments in the focus point 
scenario of the minimal SUGRA model corresponding to large sfermion masses and 
moderate to large tan f3. There is a viable region of moderate fine-tuning in the pa- 
rameter space, around tan (5 ~ 5, where the experimental limits on these electric dipole 
moments can be satisfied without assuming unnaturally small phase angles. But the 
fine-tuning constraints become more severe for tan (3 > 10. 
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It has been long recognised that the experimental limits on the electron and neutron 
electric dipole moments (EDM) imply stringent constraints on the minimal supersymmetric 
standard model (MSSM) and in particular the minimal supergravity (SUGRA) model PJ. In 
order to satisfy these limits one has to assume either unnaturally small phase angles (< 10~ 2 ) 
in the model or multi-TeV superparticle masses ||. More recently it has been shown by the 
authors of ref. || that the problem is alleviated to a large extent by internal cancellation 
between different supersymmetric (SUSY) contributions to these EDMs. Consequently, one 
can satisfy the experimental constraints on the EDMs with moderate phase angles and 
moderate superparticle masses in the unconstrained version of the MSSM |4], ||. However 
one still requires large superparticle masses in the minimal SUGRA model || §|, 0], which 
is undesirable for three reasons. It implies, i) a large fine-tuning parameter for radiative 
breaking of electroweak symmetry, ii) a less viable SUSY signal at the forthcoming colliders, 
and iii) a very large dark matter density of the universe J7| . 

Recently, Feng, Matchev and Moroi || have pointed out that the radiative electroweak 
symmetry breaking condition and hence the resulting fine-tuning is practically independent 
of the universal soft scalar mass parameter mo in the minimal SUGRA model for moderate 
to large values of tan/3( ^ 5). This is also the range favoured by the LEP data[[J. This 
has been referred to as the focus point phenomenon. It implies that one can have a m and 
hence sfermion masses of the first two generations in the multi-TeV region without affecting 
the fine-tuning parameter of electroweak symmetry breaking. Besides, one expects in this 
case an inverted hierarchy of squark masses, resulting in a distinctive SUSY signal at the 
LHC from gluino production llOL [T]J . Moreover, it has been shown to predict a dark matter 



density of the universe, which is in fact in the desired range (12 



In this paper we have calculated the electron and neutron EDMs in the focus point 
scenario to see if they can be reconciled with the corresponding experimental limits without 
assuming unnaturally small phase angles. The large mass of the first generation sfermions 
in this model helps to suppress the electron and neutron EDMs. Moreover, a large value 
of the trilinear coupling parameter Aq helps to suppress them further via a more effective 
cancellation between the different SUSY contributions. But this is partly offset by the 
increase of these EDMs with tan (3. Thus for tan j3 > 10, one cannot satisfy the experimental 
limits without assuming unnaturally small phase angles. However, there is a viable region 
of the parameter space at around tan/3 ~ 5, where the experimental limits can be satisfied 
with moderate values of the phase angles. 

In the following section we briefly discuss the focus point scenario of the minimal SUGRA 
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model and estimate the fine-tuning parameter over the region of interest to the EDM calcula- 
tion. In the next section we discuss the EDM calculation and identify the region of parameter 
space, where the EDMs can be reconciled with the experimental limits for moderate phase 
angles. We shall conclude with a brief summary of our results. 

Focus Point and Fine-tuning: 

The basic parameters of the minimal SUGRA model are mo, M1/2, Aq, B and \x - i.e. soft 
supersymmetry breaking scalar and gaugino masses, trilinear and bilinear couplings, along 



with the supersymmetric Higgs mass parameter |13fl . The last two can be determined in 
terms of the two Higgs vacuum expectation values, v\ and t>2, using the two minimisation 
conditions. The first condition determines B in terms of the ratio 1)2/ v\ = tan/3 and the 
sum 

v 2 = v \ + v 2 = 2m 2 z /(g 2 + g' 2 ) ~ 175 GeV. (1) 
The second condition gives 

l -m% = m *-f* tan2P -\rf + A R , (2) 



where the last term comes from the radiative correction to the Higgs potential. This equation 
determines the modulus of fi. 

Thus for any tan (3, the naturalness of the electroweak symmetry breaking scale requires 
Tn 2 H2 and | /j, | 2 to be of the order of m 2 z , so that there is no large cancellation between these 
quantities flUf . Since m^ 2 is linearly related to m^M 2 ^ and \A \ 2 y i a hs renormalisation 



group equation (RGE), one usually assumes the naturalness criterion to imply each of these 
parameters to be < 1 TeV. Indeed, most of the phenomenological studies within the minimal 
SUGRA model are based on this assumption. However, as pointed out by Feng et al J5|, for 
physical values of the top quark mass and the gauge couplings, m 2 ^ at the electroweak scale 
becomes practically independent of its GUT scale value tUq for tan (3 ^ 5. One can see this 
result from the approximate analytic solution of the one- loop RGE for rn 2 H [15, |TjJ. For 
tan j3 not too large, one gets while neglecting the b Yukawa coupling contribution, 

3 

m H 2 - m l ~ 2 ym o + /( M i/2> A), y), (3) 

Here / is a quadratic function of the soft parameters M1/2 and Aq, and y represents the top 
Yukawa coupling squared relative to its fixed point value, i.e. 

h\ 1 + 1/ tan 2 /? 
V ~ h 2 ~ l + l/tan 2 /3/ ^ 
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The top Yukawa coupling is related to its running mass, 



ht = m t (M t )/v sin /3, 



(5) 



which is related in turn to the physical top quark mass via 



M t = m t (M t ) [1 + Aqcd + Asusy] 



(6) 



The QCD and SUSY radiative corrections add about 6% and 4% respectively to the running 
mass to arrive at the physical top pole mass, M t = 175 ± 5 GeV 0. It is well known now 
that a physical top mass of 175 GeV corresponds to the fixed point value, tan/5/ ~ 1.5 at 



the electroweak scale |16fl , which defines the minimal value of tan/3 in this model. Such a 
low value of tan f3 is of course ruled out by the recent LEP limit on the lightest higgs boson 
mass ||, suggesting tan/3 > 2(4) for maximal (small) stop mixing. Substituting the above 
value of tan/?/ in (|) gives 

y~ 2/3 for tan/?>5. (7) 
Thus over a large range of tan (3, which is also favoured by the above mentioned LEP data, 



ni 



\ 2 of eq. (H) at the electroweak scale is practically independent of its GUT scale value 



m 



This is the so called focus point phenomenon, which implies that ttlq can be made > 1 TeV 
without affecting the naturalness criterion. The corresponding squark and slepton masses of 
the first two generations remain large at the electroweak scale, 



m 2 - T ^ 

q,l 



m 2 + 0{Ml 2 ) > 1 TeV 



since their RGEs are not affected by the top Yukawa coupling. Interestingly, the focus point 

2 

condition ensures that at the electroweak scale is practically the same as its GUT scale 



value, since 16 



|//| 2 = l.8|/i | 2 (l-y)5 



(9) 



The sensitivity of the electroweak scale to the SUSY parameters are determined from eq. 
in terms of the partial derivatives 



a, 



a dm 2 7 



m| da 



(10) 



where a denotes mo, Mi/ 2 , yUo and A . The fine-tuning is defined by the largest of these 
quantities |14 

C = max |C mo , Cm 1/2 , C^, Ca } • (11) 
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This parameter is a plausible though not unique measure of fine-tuning. It is based on the 
sensitivity of the electroweak scale to the SUSY parameters, but not other quantities like 

m t @. 

For estimating the fine-tuning parameter we have taken the radiative electroweak sym- 
metry breaking code of ref. JT^], which uses two-loop RGEs along with two- loop QCD 
correction to the top quark mass of eq. (|6|); and added the one-loop SUSY correction to the 
latter following ref. [JI8] . The radiative correction to the Higgs potential in (|2|) is evaluated 
using the complete one-loop result |19| . 



We have computed the fine-tuning parameter C in the (m — A ) planes of Figs, (la) to 
(Id) for fixed values of M\i 2 . Figs, (la) and (lb) show the contour plots of C for Mi/ 2 = 300 
and 500 GeV at tan (3 = 5, while Figs, (lc) and (Id) show the analogous plots at tan (3 = 10. 
The phases in these figures correspond to = 0, and 4>a = and 7r, for the upper and the 
lower half regions of the contours respectively. In general, C is very weakly sensitive to the 
phase cf)^ as long as M is in a range which satisfies the EDM constraints in a broad region 
of parameter space. On the other hand, it is modestly sensitive to the phase <$a as can be 
seen by comparing the 4>a = and n parts of the contours in Figs, (la) to (Id). 

We see that for moderate values of \A \ ( < 2000 GeV ) the fine-tuning parameter 
increases appreciably with M1/2, but it is effectively independent of tuq at fixed M\i 2 and 
Aq. Figs, (la) and (lb) indicate that, for contours with <pA = 0, one can go up from 
Aq — mo ~ to mo — 2000 GeV and Aq ~ 1500 GeV without paying any appreciable 
price in terms of fine-tuning. Finally, these figures also show that there is only a marginal 
improvement of the fine-tuning parameter in increasing tan/3 from 5 to 10. 

EDMs of Electron and Neutron: 

The EDM of an elementary fermion (electron or quark) is the coefficient d$ of the effective 
Lagrangian 

C E = ^d^a^F^, (12) 

which has the nonrelativistic limit dfi/j A a ■ Ei/ja, being the large component of the Dirac 
field. Fig. (2) shows the one-loop contributions to the effective Lagrangian of eq. (|12D in 
the MSSM, coming from the chargino and the neutralino exchanges, along with the gluino 
exchange in the case of quark. Denoting the generic interaction Lagrangian by 

- C int =J2M K ikPL + L ik P R )^4 h + H.C, (13) 

ik 
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the one-loop EDM is given by 0] 



It ( 47r ) m k 



QiA ( 41 + Q*s (4 



(14) 



where Pr,^ = (1 =F 7s)/2 and 



1 /o 2 ^ r \ r,/ x 1 A 2r£nr\ 

Here Q denotes electric charge. The Q t and terms in (|14D correspond to the diagrams 
with photon coupling to the chargino xt and the sfermion respectively. 

The presence of CP violating phases in the minimal SUGRA model is responsible for 
a nonzero imaginary part for the product K^L\ k in ([14]). If one neglects sfermion flavour 
mixing to avoid large flavour changing neutral current effects, then there are two independent 
physical CP violating phases in this model [p0[| . They can be chosen to be the phases of // 
and A , namely </> M and 4>a , while M 1 / 2 and Bfi are chosen to be real P, KL H|. The reality 
of B\i ensures that the Higgs vacuum expectation values and the resulting tan/5 are real. 
Following the renormalisation group equation of /i one may note that the phase of fi is scale 
independent. 

The form of the effective Lagrangian of eq. flT^D requires different chiralities of the initial 
and the final state fermions, as indicated in Fig. (2) and eq. fllil) . For the gluino exchange 
contribution, this comes from the chirality flip of the sfermion via the L-R mixing term in 
its mass squared matrix. For the chargino exchange contribution, this is accomplished via 
gaugino-higgsino mixing in the x ± mass matrix, while the sfermion preserves its chirality. 
The neutralino exchange receives contribution from both of these sources. Since both the 
L-R mixing sfermion mass squared term and the higgsino-sfermion-fermion coupling are 
proportional to mj, all the contributions are proportional to the external fermion mass. 
Another consequence of the chirality flip is the explicit proportionality of the contributions 
to the exchanged fermion mass in eq. (|H|). 

The gluino exchange contribution to the quark EDM is given by 

k=i qk \ ih / 

where D q is the L-R mixing matrix for the squark q, which diagonalises the corresponding 



* * Tfl 

Im(D q 21 D q u ) = -Im(Dl 2 Dl 2 ) = 9 {\A q \ sin0, + |p| sin0 M ig, (17) 

M q 1 ~ M q 2 
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where 

R u = cot/?, R d = tan/3, (18) 

and <j) q is the phase of A q at the electroweak scale. For the first generation of fermions, the 
magnitudes and phases of the A parameters at the electroweak scale are close to those of Aq 
at large A , since [f| 

A u ~ (1 - 0.hy)A o - 2.8M 1/2 
A d ~ A -3.6M 1/2 
A e ~ A -0.7M 1/2 . 



The chargino exchange contributions to the EDMs are given by 



—eam~ 



A\f2-Kmw sin 2 6w sin /? 



1=1 



-earrid 



Ay/2- 



wm w sin 6> w cos {3 i=1 

2 



EMM 




il- 



eum, 



m 



m 



A\/2m w sin 2 6 W cos /? i=1 



Aft 



(19) 

(20) 
(21) 
(22) 



where Wi,<ii refer to the dominantly left-handed squark mass eigenstates. In our numerical 
analysis we have also included the small contributions from the terms with w 2 , <i 2 - Here, 
U and V are the gaugino-higgsino mixing matrices, which diagonalise the chargino mass 
matrix. Explicit expression for the U and V matrices are given in ref. in terms of Mi/ 2 , 
tan/3, \fjj\ and M . We shall simply note here that each of the coefficients Im(U* 2 V/i) and 
Im{y^U*]), is proportional to |/i| sin0 M . Consequently, 



d 



u,d,e 



oc \y,\ sin^ 

The neutralino exchange contributions to the EDMs can be collectively expressed as 

ea 



(23) 



a xo 



An sin 2 9 



2 4 m 



W k= i i=i 



fk 




(24) 



where 



/ 
Vik 



V2{tanM<2/ " T 3/)^ + T 3/iV 2i } D(l - K/iV w D£ 



f 

lk I 5 



(25) 
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with 

m u TTl^e 
= /— 5 e = /— V J 

V2mvysin/3 v2mjycos/3 
and b = 4(3) for u(d,e). The are the L-i? mixing matrices for the sfermion /, which 
occurred earlier in eq. (0). Explicit expression for its matrix elements are given in ref. || 
in terms of \Af\, |//|, 0/ and M . Finally, N is the 4x4 unitary matrix, diagonalising the 
neutralino mass matrix, which is evaluated numerically. 

The main contributions to the EDM of quarks come from chargino and gluino exchanges, 
while neutralino exchange contribution is relatively small. They are related to the neutron 
EDM via the nonrelativistic quark model relation |2~lH , 

1 



d r 

3 



Ad d - d u 



Ve, (27) 



where tje = 1-53 is a QCD correction factor for evolving down the quark EDMs from the 
electroweak to the hadronic scale || 

There are two other contributions to the neutron EDM, arising from the quark chromo- 
electric dipole moment and the gluonic dimension-six operator, which are defined by the 
effective Lagrangians 

C'b = ^dlqa, ul ,T a qG^ a (28) 

and 

£-E = ~ Q^C] C abcG fivaGl' °G^ C , (29) 

where T a are the SU(3) generators, f a bc the Gell-Mann coefficients and G^ ua the gluonic field 



tensors [23[] . Their contributions to the neutron EDM were earlier supposed to be small with 
respect to the quark EDM contribution of eq. (|27| ) |22|j . But as demonstrated in ref. [[J , the 
large internal cancellation between the chargino and the gluino contributions to the d q can 
make the net quark EDM contribution comparable to those from the quark chromoelectric 
and the gluonic dimension-six operators over certain regions of parameter space. In the 
present analysis we have included each of these three contributions to the neutron EDM 
following ref. 

We have investigated the ranges of the phase angles <fi^ and <pA Q , over which the predicted 
electron and neutron EDMs can be reconciled with the corresponding experimental limits 
I, 

d e < 4.3 x KT 27 ecm, d n < 6.3 x KT 26 ecm. (30) 

It may be noted here that the above neutron EDM limit is a factor of 2 smaller than that 
considered in most of the previous analyses || || 0. Although it is still an order of 



S 



magnitude larger than the electron EDM, both of these provide comparable constraints over 
the parameter space of interest in our analysis. 

Figs. (3a), and (3b) show the range of at tan/3 = 5, over which the predicted 
electron and neutron EDMs can be reconciled with the experimental limits fl30|), by allowing 
a variation of 4>a over the range of — 7r to tt. It should be mentioned here that the chargino 
contribution to d n or d e overshoots the corresponding experimental limits for the moderate 
values of |0 M | (~ 0.1 — 0.2 radian) shown in Figs. (3a) and (3b). What helps to satisfy the 
experimental limits is a cancelling contribution from gluino exchange for d n (and neutralino 
exchange for d e ). Consequently, there is a strong correlation between the two phase angles, 
as noted in earlier analyses. In particular the maximal allowed value of for a given m 
and \Aq\ corresponds to \4>a \ ~ tt/2, and there is an opposite sign correlation between the 
phases. The relatively smaller range of M at moderate \Aq\ (< 2000 GeV) is due to a larger 
coefficient of sin M in the chargino contribution of eq. (|23| ) in comparison with the coefficient 
of sin^Q in the gluino contribution coming from eqs. (|l6|JT7|jr9|) . It is seen that only for 
very large \Aq\ (> 6000 GeV) the two coefficients become comparable; and one can satisfy 
the experimental limits for any value of M . But, one has to pay a high price in terms of the 
fine-tuning parameter amounting to C > 1000. Besides, the purely gluonic dimension six 
operators play an effective role in the cancellation mechanism in this region. 

It should be added here that in plotting Figs. (3a) to (3d) we have scanned the mo space 
in 50 GeV bins. Decreasing the size of this bin further leads to occasional spikes in the 
maximum value of 10^]. This is indeed an important effect reflecting further suppression of 
the neutron EDM due to internal cancellation amongst the electric dipole, chromoelectric 
dipole and the gluonic operator contributions. As already noted in ref. ||, such a drastic 
suppression of the EDM can occur over narrow ranges of the SUGRA parameters due to 
these cancellations. But in the present analysis we shall concentrate only on those results 
which hold over wide range of SUGRA parameters, granting a possible correlation between 
the two phases. 

In Figs. (3a) and (3b) we have indicated the fine-tuning parameter C at some specific 
points on the fixed \A Q \ contours. Essentially, the point marked on each contour corresponds 
to the region of the right tip of the fixed C contours for <pA — of Figs, (la) and (lb). 
Comparing the two figures one can easily see that C hardly varies up to the marked point on 
each \Aq\ contour. Thus one can accommodate at least moderate values of l^™ 5 ™] = 0.1 — 0.2 
radian for C = 100 — 200, i.e. without paying any price in term of fine-tuning. Figs. (3c) 
and (3d) show the analogous ranges of M at tan/? = 10. We see that for given values of 
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\Aq\, corresponding to the similar values of C, the M range is less than half the size of that 
of Figs. (3a) or (3b). The underlying reason of course is the comparatively larger coefficients 
of the chargino contributions for d d of eq. ( ETj ) and d e of eq. (E2I) at large tan/3. Thus the 
EDM limits disfavour large values of tan/? (> 10). 

Finally, we also analyse the case of tan/? = 3 and My 2 = 300 GeV, as displayed in 
Figs. (4a) and (4b). Fig. (4a) shows the contours for constant fine-tuning C which are very 
different from what we found in Fig. (la). Such a low value of tan/? falls outside the focus 
point scenario ||. Besides, it is disfavoured by the LEP limit on the lightest Higgs mass [f|. 
Nonetheless most of the previous EDM analyses have concentrated in this low tan /? region 
H H |?J , since it corresponds to smaller coefficients of the chargino contributions of eqs. fl2~l~| ) 
and (|22|). However, in this region of tan/? the fine-tuning parameter C steadily increases 
with mo unlike what one finds in the focus point scenario. Here, C is same as C^ ; and a 



I) I 'l i' 1 1 lit \ i ir 1 1 ' 



contour of constant C^ is a part of an ellipse [17] • Fig. (4b) shows the variation of maximal 



with m for various \A \ values. Unlike Figs. (3a) to (3d), C increases here rapidly 
along the contours of constant \Aq\. Consequently, a |0™ ax | of 0.1 radian would correspond 
to a fine-tuning measure C ~ 200, which is larger than the value required at tan/? = 5 (Fig. 
(3a)). 

Summary: 

We have analysed the electron and neutron EDMs in the focus point scenario of the min- 
imal SUGRA model along with the fine-tuning parameter. In this scenario the soft scalar 
mass mo can go up to 2 TeV without affecting the fine-tuning parameter. Similarly, the 
trilinear coupling parameter can be increased from to 1.5 TeV without any appreciable 
increase in fine-tuning. The large m values correspond to large masses for the 1st generation 
of sfermions which helps to suppress the one-loop SUSY contributions to the EDMs. More- 
over, the large |^4o| corresponds to larger gluino (neutralino) contribution to quark (electron) 
EDM, which can cancel the chargino contribution more effectively. Therefore, one can satisfy 
the experimental limits of the electron and the neutron EDMs without assuming unnaturally 
small phases (j)^ and 4>a for mo and \A \ values of ~ 2 TeV each. But this is possible only 
for a moderate value of tan/? ~ 5. The chargino contributions to the EDMs increase with 
tan/?, so that the experimental limits cannot be satisfied without assuming small \<p^\ or a 
large fine-tuning parameter C for tan/? > 10. Since the completion of this work a general 
phenomenological analysis in the focus point scenario including the EDMs has appeared 
recently in ref. |[24|| . However the present work contains a more detailed treatment of this 
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issue. 

U.C. and D.P.R. thank Probir Roy for helpful discussions. 
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Figure 1: Lines of constant fine-tuning C in (mo — Aq) plane for different values of tan f3 and 
Mi/2- Here 0^ = 0, and for each contours 0a o corresponds to and 7r for the upper and the 
lower parts respectively. The shaded areas in the right represent the excluded regions due 
to the chargino mass limit and the electroweak radiative breaking constraint. The shaded 
areas in the left are excluded by the top squark mass bounds. 
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Figure 2: One loop diagrams contributing to the electric dipole moments of quarks and 
leptons. 
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Figure 3: Lines of constant \A \ in the plane of (\4>^ max \ — m ) for different cases of tan/5 
and Mi/ 2 . Here the maximum value of is obtained by varying <f>^ from — n to n. The 
symbols shown in the figures refer to the appropriate values of the fine-tuning measure C 
which remains practically constant on each curve up to the marked point. 
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Figure 4: (a): Lines of constant fine-tuning C in (m — A ) plane for tan/5 = 3 and 
M1/2 — 300 GeV. Here = 0, and for each contours 4>a corresponds to and tt for the 
upper and the lower parts respectively. The shaded areas in the left are excluded by the top 
squark mass bounds, (b): Lines of constant \Aq\ in the plane of \(f) l / nax \ — mo for tan/5 = 3 
and Mi/2 — 300 GeV. Here the maximum value of is obtained by varying <pA from —n 
to 7T. Unlike Figs. (3a) to (3d), here C strongly varies along the constant \A \ contours. 
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